Mesoscopic to universal crossover of transmission phase of multi-level quantum dots 
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Transmission phase a measurements of many-electron quantum dots (small mean level spacing 8) 
revealed universal phase lapses by tv between consecutive resonances. In contrast, for dots with only 
a few electrons (large 8), the appearance or not of a phase lapse depends on the dot parameters. 
We show that a model of a multi-level quantum dot with local Coulomb interactions and arbitrary 
level-lead couplings reproduces the generic features of the observed behavior. The universal behavior 
of a for small 8 follows from Fano-type antiresonances of the renormalized single-particle levels. 

PACS numbers: 73.23.-b, 73.63.Kv, 73.23.Hk 



One of the longest-standing puzzles in mesoscopic 
physics is the intriguing phase-lapse behavior observed in 
a series of experiments [H S 13 on Aharonov-Bohm rings 
containing a quantum dot in one arm. Under suitable 
conditions in linear response, both the phase and magni- 
tude of the transmission amplitude T = \T\e ta of the dot 
can be extracted from the Aharonov-Bohm oscillations 
of the current through the ring. If this is done as func- 
tion of a plunger gate voltage V g that linearly shifts the 
dot's single-particle energy levels downward, Ej = e° — V g 
(j = 1,2,... is a level index), a series of well-separated 
transmission resonances [peaks in |T(V^)|, to be called 
"Coulomb blockade" (CB) peaks] of rather similar width 
and height was observed, across which a(V g ) continu- 
ously increased by n, as expected for Breit-Wigner-like 
resonances. In each CB valley between any two succes- 
sive CB peaks, a always jumped sharply downward by tt 
("phase lapse", PL). The PL behavior was observed to 
be "universal" , occurring in a large succession of valleys 
for every many-electron dot studied in [l|, 2, 3]. This uni- 
versality is puzzling, since naively the behavior of oe(V g ) 
is expected to be "mesoscopic", i.e. to show a PL in 
some CB valleys and none in others, depending on the 
dot's shape, the parity of its orbital wavefunctions, etc. 
Despite a large amount of theoretical work (reviewed in 
[US]), no fully satisfactory framework for understanding 
the universality of the PL behavior has been found yet. 

A hint at the resolution of this puzzle is provided by 
the most recent experiment Q, which also probed the 
few-electron regime: as V g was increased to successively 
fill up the dot with electrons, starting from electron num- 
ber iV e = 0, a(V g ) was observed to behave mesoscopically 
in the /ew-electron regime, whereas the above-mentioned 
universal PL behavior emerged only in the many-electron 
regime (N c > 15). Now, one generic difference be- 
tween few- and many-electron dots is that the latter have 



Tj stemming from hybridization with the leads. Thus, 
Ref. l3j suggested that a key element for understanding the 
universal PL behavior might be that several overlapping 
single-particle levels simultaneously contribute to trans- 
port. Due to the dot's Coulomb charging energy U, the 
transmission peaks remain well separated nevertheless. 

Previous works have studied the transmission ampli- 
tude of multi-level, interacting dots ESBIIGli How- 
ever, no systematic study has yet been performed of the 
interplay of level spacing, level widths, and charging en- 
ergy that combines a wide range of parameter choices 
with an accurate treatment of the correlation effects in- 
duced by the Coulomb interaction. The present Let- 
ter aims to fill this ga p by using two powerful methods, 



the numerical (NRG) [llj, [lj and functional (fRG) [R 



renormalization group approaches, to study systems with 
up to 4 levels (for spinless electrons; see below). We find 
that if the ratio of average level spacing 6 to average 
level width T is decreased into the regime 6 < T, one 
of the renormalized effective single-particle levels gener- 
ically becomes wider than all others, and hovers in the 
vicinity of the chemical potential fj, in the regime of V g 
for which the PLs occur. Upon varying V g , the narrow 
levels cross /i and the broad level, leading to Fano-type 
antiresonances accompanied by universal PLs. For 8 > T, 
ct(V g ) behaves mesoscopically 1J] for all U. Decreasing 
5 thus causes the PL behavior to generically change from 
mesoscopic to universal, as observed experimentally 
Model: The dot part of our model Hamiltonian is 
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3+3 



with 



= d)dj 



and dot creation operators d] for spin 



less electrons, where U > describes Coulomb repulsion. 
The semi-infinite leads are modeled by a tight-binding 



-j+i 



smaller level spacings in- 
filled levels. With increasing jV e , their <5j's should even- 
tually become smaller than the respective level widths 



e° for the topmost chain H i = -*Em=o( c L/">+i,l + H.c.) and the level- 



lead couplings by Ht = —J2ji.(t l j c oidj + H.c), where 
Cm j annihilates an electron on site m of lead I — L, R and 
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tj are real level-lead hopping matrix elements. Their rel- 
ative signs for successive levels, Sj = sgn(tft^tj +1 t^ +1 ), 
are sample-dependent random variables determined by 
the parity of the dot's orbital wave functions. The ef- 
fective width of level j is given by Tj = Tf + Tf, with 
= irp\t l j\ 2 . We take p, the local density of states at 
the end of the leads, to be energy independent, choose 
p = 0, and specify our choices of tj using the notation 

a = {si,s 2 , . . .}, 7 = {rf , rf , r£, . . .}, r = £ jV r< . 

Methods: We focus on linear response transport and 
unless stated otherwise, on zero temperature (r = 0). 
Then the dot produces purely elastic, potential scattering 
between left and right lead, characterized by the trans- 
mission matrix T w = 2-wp £) y t\Qfj (0)t l j , where Gfj(w) 
is the retarded local Green function which we compute 
using NRG and fRG. The NRG is a numerically exact 
method that is known to produce very accurate results 
Hill. The fRG s a renormalization procedure for the 
self-energy £ and higher order vertex functions (see (l3j 
for details). We use a truncation scheme that keeps the 
flow equations for S and for the frequency independent 
part of the effective two-particle (Coulomb) interaction. 
Comparisons with NRG [l3( have shown this approxima- 
tion to be reliable provided that the number of (almost) 
degenerate levels and the interaction do not become too 
large. fRG is much cheaper computationally than NRG, 
enabling us to efficiently explore the vast parameter space 
relevant for multi- level dot models. 

At the end of the fRG flow, the full Green func- 
tion at zero frequency takes the form [£7^(0)]^ = 
— hij + iAij, with an effective, noninteracting (but V g 
and [/-dependent) single-particle Hamiltonian hij — 
(e® — V g )Sij + Sjj, whose level widths are governed by 
TrpJ2i^ l j- To interpret our results, we adopt the 
7. , with eigenvalues — Sj + iTj, and 
j as level positions and widths of a renor- 
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and r 



eigenbasis of [^' R (0)] i 



view e 

malizcd effective model (REM) describing the system 

For LR-symmetry, T^ — T^, an NRG-shortcut can be 
used, which is much less demanding than computing the 
full Gfj{u): the S-matrix is then diagonal in the even-odd 
basis of the leads and its eigenvalues depend on the total 
occupancies n± of all levels coupled to the even/odd lead 
(Friedel sum rule), so that the transmission amplitude 
T = T LR takes the form T = sin[7r(n+ - n _)] e " r (™++™-). 
A transmission zero (TZ) and hence PL occurs when 
n + = n_ modi [Figs. [lje,h): n± in thin dashed/dashed- 
dotted line]. 

Results: Our results are illustrated in Figs. [1] to [31 
fRG and NRG data generally coincide rather well (com- 
pare black and orange lines in Figs. [T] and [2]), except 
for N — 4 when bothC/^r,<5<r, and correlations 
become very strong [Fig. EJf )] . The figures show the fol- 
lowing striking qualitative features, that we found to be 
generic by running the fRG for ten-thousands of parame- 
ter sets, which is possible as a complete T(V g ) curve can 



13 



N_2 x=0 y a= /r={0.3,0.3,0.7,0.7} y n=+ /r={0.2,0.4,0.4,1} 



u/r=o.4 



(a) 



u/r=6 



(b)~- n NRG (y ) ^-="(C) 

.... n b , NRG (t o _ J / 



u/r=30 



in 



j 5/1^10: ME SOS CO PIC \ 

A I i is I I 




FIG. 1: |T(V S )| and a(V g ) for N — 2, s§,i = ±5/2, and r = 0: 
decreasing S/T produces a change from (a,b,c) mesoscopic via 
(d,e,f) crossover to (g,h,i) universal behavior; increasing U/T 
leads to increased transmission peak spacing. (b,e,h) include 
the occupancies n+ (thin dashed) and n_ (thin dash-dotted) 
of the levels coupled to the even/odd lead in the case of LR- 
symmetry. The condition n+ = n_ mod 1 produces a TZ and 
PL. For the blip and hidden TZ near V g = in (i), see pj . 

be obtained within a few minutes on a standard PC: 

Mesoscopic regime: For S > T [Figs.[TJa,b,c),[5]Ja,b,c)], 
we recover behavior that is similar to the U — case. 
Within the REM it can be understood as transport oc- 
curring through only one effective level at a time [see 
Fig. [3][a,b,c)], with Tj ~ Tj. Each ij that crosses p 
produces a Breit-Wigner-like transmission resonance of 
width 2Tj and height governed by T^ /T^. At the cross- 
ing the other levels are shifted upwards by U [charging 
effect; Fig. [3ta)] leading to renormalized peak separa- 
tions ("level spacings") Sj + U . Between two peaks, 
a(Vg) behaves mesoscopically: depending on the sign Sj 
one either observes a PL (sj = +) or continuous evolu- 
tion of a (sj = — ) (l4j |. Additional PLs occur to the left 
or right, beyond the last transmission peak [Fig.^a)]. 

Mesoscopic to universal crossover: As the ratio S/T 
is reduced, the behavior changes dramatically: the TZs 
and PLs that used to be on the far outside move inward 
across CB resonances [see evolution in Figs. [IJb,e,h)]. 

Universal regime: A universal feature (ljj emerges for 
^ L < U (crossover scales are of order 1, but depend 
on the chosen parameters): for all choices of the signs a 
and generic couplings 7, the N CB peaks over which a 
increases by tt are separated by N — 1 PLs, each accom- 
panied by a TZ [Figs. QJg,h,i), [2{d,e,f)]. This is consis- 
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N=4 



U/r=2, a={+-+}, 



y/r={.48,1.12,.4,1.2,.16,.24,.12,.28} 



u/r=4, a={— }, x/r= 

y/r=(.266,.533,.4,.8,.133,.266,.533,1.07} 



FIG. 2: \T(V a ; 




U/T=24, a={--+}, T/r=0.16 
yr={.16,.24,.32,.48,.72,.48,.68,.92} 



for N — 4, with equidistant levels 



The qualitative features do not change if this assumption 
is relaxed, or if U is assumed to be slightly level-dependent, U — * Ujji. Decreasing S/T produces a crossover from (a,b,c) 
mesoscopic to (d,e,f) universal behavior; increasing U/T leads to increased spacing of the transmission peaks and PLs. For 
clarity the finite temperature (r > 0) curves were shifted downwards (by 0.1 for |T|). (f) For U/T 3> 1 and t/T — 0.16, the 
CB peak and PL shapes are strikingly similar to those observed experimentally at comparable ratios of t/T (see Fig. 3 of 0] 
and Fig. 6 of [3j). For small 8/T and U/T ^> 1 fRG becomes less reliable and the results begin to differ from those of NRG. 



tent with the experimentally observed trend. For small 
to intermediate values of U/T [Figs. DIg,h), I21d,e)], the 
transmission peaks are not well-separated, and o:(V g ) has 
a saw-tooth shape. As U /T increases so does the peak 
separation and the corresponding phase rises take a more 
S- like form [Figs. QIh,i) , e,f )1 . At finite temperatures of 
order r > 5 [17[ sharp features are smeared out (Fig. \2§ ■ 
For U /T as large as in Figs.QJi) and[^f), the behavior of 
a(V g ) (both the S'-like rises and the universal occurrence 
of PLs in each valley) as well as the one of |T(V S )| (simi- 
lar width and height of all CB peaks) is very reminiscent 
of that observed experimentally, in particular for r > S 
[Fig. Of)]. For U/T » 1 the full width of the CB peaks 
is of order 2 AT (not 2Tj as in the mesoscopic regime), 
indicating that several bare single-particle levels simul- 
taneously contribute to transport. The r dependence 
of the width of the PLs is different from the behavior 
t 2 / (5 + U) 2 found in the mesoscopic regime [l^ and will 
be discussed in an upcomming publication. Note that for 
the temperatures considered here the width of the PLs is 
still much smaller than the width of the CB peaks. 

For certain fine-tuned parameters (7 and a) the behav- 
ior at small <5/r deviates from the generic case. For = 
2 the nongeneric cases were classified in [9]. In Fig. [l]only 
generic parameters are shown. For N > 3 LR-symmetric 
couplings produce nongeneric features. However, these 
features are irrelevant to the experiments. They quickly 
disappear upon switching on LR-asymmetry or r > 0. 

Interpretation: We can gain deeper insight into the ap- 
pearance of the TZs and PLs in the universal regime from 
the properties of the REM obtained by fRG for moderate 
U/T [at which NRG and fRG agree well; Fig.^e)]. For 



N > 3, 5 < r and U — 0, two of the effective levels are 
much wider than the others, since A^, being a matrix 
of rank 2, has only two nonzero eigenvalues [18| (for the 
N = 2 case, see 
at U > 0: for N 



We found that this also holds 
3,4 one effective level is typically 
a factor of 2 to 3 wider than the second widest, while 
the remaining 1 or 2 levels are very narrow [Fig. |3ff)]. 
At 5 < T [Fig. EJd)] the interaction leads to a highly 
nonmonotonic dependence of Ej on V g which is essen- 
tial for our universal PL scenario: As V g is swept, the 
widest level hovers in the vicinity of \x over an extended 
range of V g values, whereas the narrow ones cross /i - and 
therefore also the widest one - rather rapidly. This leads 
to a Fano-type effect [13, [H [H, [H, [2J] whose effective 
Fano parameter q is real, by time-reversal symmetry [2l| . 
Thus, each TZ, and hence PL, can be understood as a 
Fano-type antiresonance arising (irrespective of the signs 
of t{) from destructive interference between transmission 
through a wide and a narrow level. The crossings of the 
narrow levels and /i, and thus the PLs, are separated 
by U due to charging effects. In contrast, for U = 0, 
ij oc — V g for all renormalized levels and no levels cross 
each other. Our fRG studies indicate that for the regime 
8 < F, the Fano-antiresonance mechanism is generic for 
U ^ r. We thus expect it to apply also for interactions 
U ^> r for which fRG is no longer reliable. 

The fact that the combination of a wide and several 
narrow levels leads to PLs was first emphasized in Q 
(without reference to Fano physics). However, whereas in 
[8| a bare wide level was introduced as a model assump- 
tion (backed by numerical simulations for noninteract- 
ing dots of order 100 levels), in our case a renormalized 
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H.2.1, BSF and the Humboldt foundation. 



FIG. 3: Renormalized single-particle energies Sj, occupancies 
hj and level widths Pj (dashed) of the REM, and the resulting 
|T| and a (all as functions of V g ), for N — 4 and 5j = 5 at 
t = 0. The parameters are the same as in Fig. [2jb,e). 

wide level is generated for generic couplings if d < T. 
Also, whereas in [f| the wide level repeatedly empties 
into narrow ones as V g is swept (because <C U was 

assumed), this strong occupation inversion [25j | is not re- 
quired in our scenario. In Fig. (3^e), e.g., the wide level 
remains roughly half-occupied for a large range of V g , 
but TZs and PLs occur nevertheless. We thus view oc- 
cupation inversion, if it occurs, as a side effect, instead 
of being the cause of PLs [2^| . 

Conclusions: The most striking feature of our results, 
based on exhaustive scans through parameter space for 
N = 2,3,4, is that for any given, generic choice of cou- 
lings (7 and a), the experimentally observed crossover 
from mesoscopic to universal a(V" s )-behavior can be 
achieved within our model by simply changing the ratio 
8/T from > 1 to < 1, provided that U > T. The univer- 
sal 7r PLs result from Fano-type antiresonances of effec- 
tive, renormalized levels, which arise because interactions 
cause a broad level (occuring, if 6 < T, already for U = 0) 
to be repeatedly crossed by narrow levels. A quantitative 
description requires correlations to be treated accurately. 
We expect that the main features of this mechanism carry 
over to the case of spinful electrons, since for 5 < T spin 
correlation physics (such as the Kondo effect) does not 
play a prominent role [27| . 
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